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We investigate the statistics of Iso-height lines of (2+l)-dimensional Kardar-Parisi-Zhang model 
at different level sets around the mean height in the saturation regime. We find that the exponent 
describing the distribution of the height-cluster size behaves differently for level cuts above and 
below the mean height, while the fractal dimensions of the height-clusters and their perimeters 
remain unchanged. The winding angle statistics also confirms again the conformal invariance of 
these contour lines in the same universality class of self-avoiding random walks (SAWs). 
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The study of surfaces, their static statistical proper- 
ties, as well as growth and evolution dynamics and their 
morphological properties, has been attracting an ever in- 
creasing amount of interest over the last two decades. 
A method of characterizing surfaces is by looking at their 
iso-height contour lines and the islands that are generated 
by a cut through the surface at a certain constant height. 
Such a situation can be found in the patterns exhibited by 
topographical islands and continents. These islands and 
their coastlines have some fractal features with a fractal 
dimension related to the roughness exponent a, given by 
the structure function ([h(x) — h(x + r)] 2 ) ~ r 2a , where 
h(x) is the height of interface as a function of its position. 

Theoretical modeling of the growth processes started 
with the work by Edwards and Wilkinson (EW) [1] who 
suggested that one might describe the dynamics of the 
height fluctuations by a simple linear diffusion equation. 
Kardar, Parisi, and Zhang (KPZ) [2] realized that there 
is a relevant term proportional to the square of the height 
gradient which represents a correction for lateral growth. 
The KPZ equation is given by 

: ^M = z ,v 2 / l +^ivM 2 +^(x J t). a) 

The first term on the r.h.s describes relaxation of the in- 
terface caused by a surface tension v, and the nonlinear 
term is due to the lateral growth. The noise r\ is uncorre- 
cted Gaussian white noise in both space and time with 
zero average i.e., (?y(x, t)) = and {r](x,t)r](x.' ,t')) = 
2DS d (x-x')S(t-t'). 

In (l-fl)-dimension the roughness and growth expo- 
nents were exactly obtained, a = 1/2 and (3 = 1/3 re- 
spectively [3], while for the (2+l)-dimensional case there 
are just numerical evidence and predictions [4]. 

The KPZ equation is invariant under translations along 
both growth direction and perpendicular to it, as well as 
time translation and rotation. Despite these symmetries, 
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the existence of driving force perpendicular to the in- 
terface (due to the nonlinear term), breaks the up-down 
symmetry (h — > —h) [3]. In two dimensions these sym- 
metries and growth dynamics can affect the statistics of 
the iso-height lines (island coastlines) at different level 
sets, which is the main subject of the present work. 

In our previous paper [5] we have focused on zero iso- 
height lines of the (2+l)-dimensional KPZ model in the 
saturation regime (mean height was set to zero). Us- 
ing the theory of Schramm-Loewner evolution (SLE), 
we have numerically shown that the contour lines of 
zero height behave statistically like self avoiding walks 
(SAWs) and they can be defined by the family of confor- 
mally invariant curves i.e., SLE K curves with diffusivity 
k = 8/3. The statistics of these objects for the EW model 
has been shown to be the same as the interfaces in the 
0(2) model, and can be described by SLE4. 

The SLE process, introduced by Schramm [6] describes 
the scaling limit of a variety of statistical mechanical 
models in two dimensions (some review articles are given 
in [7]). Schramm and Sheffield showed that the contour 
lines in a two-dimensional discrete Gaussian free field are 
statistically equivalent to SLE4 [8]. Moreover, it is shown 
that the restriction property only applies in the case for 
k = 8/3 [9]. Since self-avoiding random walk (SAW) sat- 
isfies the restriction property, it is conjectured that in 
the scaling limit it falls in the SLE class with k = 8/3 
[10]. The scaling limit of SAW in the half-plane has been 
proven to exist [11] but there is no general proof of its 
existence. 

The theory of SLE has recently been applied to many 
experimental and physical systems. It is shown that the 
statistics of the zero-vorticity lines in inverse cascade 
of two dimensional (2D) Navicr-Stokcs turbulence is 
conformally invariant and belongs to the percolation 
universality class [12]. The same issue has been studied 
for zero-temperature isolincs in the inverse cascade 
of surface quasigeostrophic turbulence [13], domain 
walls of spin glasses [14] and the nodal lines of random 
wave functions [15]. Moreover, it has been shown 
recently that the statistics of the iso-height lines on 
the experimentally grown WO3 surface is the same as 
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FIG. 1: (color online). Log-log plot of the average area of a 
cluster M versus the radius of gyration R, at three different 
cuts with 5 = -2, and +2. 



domain walls statistics in the critical Ising model [16]. 
Avalanche frontiers in sandpile models have been shown 
to be conformally invariant and in the same universality 
class of loop erased random walks [17]. 

Here we briefly review the results obtained in [5] for 
zero height level cuts, and then extend them for cuts 
made at different heights. 

We have integrated the discretized KPZ equation on a 
square lattice of size 2048 2 , with periodic boundary con- 
ditions. The details of numerical integration and simula- 
tion are given in [5] . 

Consider an ensemble of 2D-KPZ saturated surfaces and 
a cut is made at specific height say hs = (h(x)) + 
Sy/([h(x) - (h(x))] 2 ) := 0, where the symbol (••) de- 
notes spatial averaging. Then we define each island 
(cluster height) as a set of connected sites with positive 
height which were identified by the Hoshen-Kopelman al- 
gorithm. 

The scaling of the mass M of a cluster with the radius 
of gyration R, behaves like M ~ R Dc , where D c is the 
fractal dimension of the cluster which is D c = 2 in this 
case. As shown in Fig. 1, this fractal dimension remains 
unchanged for different 5. 

The fractal dimension of a coastline (or loop because of 
periodic boundary conditions), can be obtained with the 
scaling relation between the average length of a loop I, 
and the radius of gyration i?, as I ~ R Dl . This fractal di- 
mension remains also fixed, within numerical errors, for 
cuts made at different 8. The best fits to data shown in 
Fig. 2, yields the fractal dimension of a contour line in 
the range of Di = 1.34 ± 0.02 (Fig. 2). 
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FIG. 2: (color online). Log-log plot of the average length of a 
loop I versus the radius of gyration R, at three different cuts 
with S = -2, and +2. 



Powerful scaling arguments made by Kondev and Henley 
[18], connect the fractal dimension of a contour line to 
the roughness exponent a of the surface, 

Di = 2-xi- a/2, (2) 

where xi is the loop correlation exponent. Although the 
exact value of x t = 1/2 is for a — and 1, [19] but it is 
conjectured that its value is super universal and is inde- 
pendent of a for Gaussian surfaces. 
In the case of 2D-KPZ surface the finite size scaling 
for the interface width yields the roughness exponent [5] 
a = 0.37 ±0.01, which is in mild conflict of Eq. (2). This 
may be because the field h(x) does not follow a Gaus- 
sian distribution. In other words, the fractal dimension 
Di = 1.34 ± 0.02 ~ 4/3 obtained for contour lines of 2D- 
KPZ surface, is equal to what one obtains for Gaussian 
surfaces with a = 1/3. 

The island size distribution has also a power-law be- 
havior. As can be seen in Fig. 3, there arc two distinct 
scaling regions for the distribution of the island size. We 
find that the small size islands are distributed according 
to a power-law distribution n(M) ~ M _Ts , with a same 
exponent t s = 2 ± 0.05 for level cuts at different 6. 

For the other region with larger island size dominant, 
the power-law behavior is held, but with different expo- 
nents at different level cuts. 

For cuts made at lower values of 6, percolative height 
clusters appear dominantly which their sizes are in the or- 
der of sample size. Inside these percolative islands there 
are some lakes (negative height clusters) which can also 
contain smaller islands with the distribution as shown in 
Fig. 3 (for 5 = —2). Increasing in the height of the cut, 
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FIG. 3: (color online). Main: log-log plot of the number 
of clusters of area between M and 1.05M, at three different 
cuts with 8 = -2, and +2. Inset: the exponents for the 
island distributions as a function of 8. The errors are less 
than 0.05 for all exponents. The slope of the dashed lines are 
— 1.77 ± 0.03, the best fitted for 8 < 0, and 0, which is drawn 
for comparison for 8 > at r = 0.75. 
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FIG. 4: (color online). Log-log plot of the number of pos- 
itive (squares) and negative (circles) height-clusters of area 
between M and 1.1M, at a level cut made at the mean height 
with 8 = 0. Solid lines are the best fits to the data in the 
linear region yielding the distribution exponent r ~ 0.77 and 
t ~ 1 for positive and negative height-clusters, respectively. 



tends to a continuous distribution of islands with differ- 
ent sizes and different scaling behavior. We find that for 
cuts below the mean height (i.e., 8 < 0), the island-size 
distribution exponent r takes different values for different 
<5, while it takes almost a same value (within statistical 
error) for positive height cuts. This behavior is shown in 
the inset of Fig. 3, which attributes two different regimes 
for level cuts i.e., 8 < and 6 > 0. For level cuts made at 
8 < 0, the island-size distribution exponent r decreases 
linearly with the slope of —1.77 ±0.03, while it takes val- 
ues around r = 0.75 within statistical errors for S > 0. 
This asymmetrical behavior of the exponent r around 
the mean height may be interpreted as the breakdown 
of the up-down symmetry under changing h — ► — h in 
the KPZ equation (1). This asymmetry in the dynamics 
of the growth process can tend to an asymmetry in the 
statistics and the distribution of valleys and overhangs. 

In a system with up-down symmetry it is expected 
that the valleys and overhangs have the same statistical 
behavior with same distribution. This behavior is not 
seen in the mean-height cuts of the 2D-KPZ surface. The 
size distribution of the islands with positive and negative 
heights differ. As shown in Fig. 4, the exponents defining 
these two distributions are quite different with values r ~ 
0.77 and r ~ 1 for positive and negative height-clusters, 
respectively. 

Although the exact cause for the behavior of the ex- 



ponent t as function of the height of the cut eludes us, 
but it may be possible to guess its behavior around the 
mean height [18] (a slightly different problem has been 
investigated in [20] for random Gaussian surfaces). 
Using scaling arguments, it is shown in [18] that the aver- 
age number density of contour lines (the coastlines here), 
scales with the radius of gyration R as n(R) ~ R~ 2+a . 
Since the dimension of the islands in our case is 2 (see 
Fig. 1), one can expect that the radius of gyration 
for the perimeter of the islands and the islands them- 
selves have the same statistical behavior. So, the aver- 
age number density of the islands of size M, is given by 
n(M) ~ M~ 1+a ' 2 . Within an uncertainty in determin- 
ing the roughness exponent a, it yields an approximated 
value for the exponent r ~ 0.81 which works here for the 
height cuts around and above the mean height. 

In the rest of the paper we investigate the confor- 
mal symmetry of the coastlines of cuts made at different 
heights of two dimensional saturated KPZ surface. The- 
ory of SLE provides an appropriate approach to check 
conformal properties of the geometrical features of such 
systems. Looking at the fractal dimension obtained for 
the KPZ coastlines at different level sets, it agrees with 
the SLE curves of fractal dimension Dj = 1 + k/8, with 
k = 8/3, conjectured to describe the scaling limit of 
SAWs. In [5], we checked various consistencies between 
the coastlines and both SAWs and SLE 8 / 3 . Such coast- 
lines can statistically be defined as the outer boundary 
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FIG. 5: (color online). Winding angle statistics of iso-height 
lines of 2D-KPZ surface simulated on square lattice of size 
2048 2 . (a) Main: logarithmic behavior of the variance of the 
winding angle as a function of distance I along the curve (see 
Eq. (3)). Upper-left inset: the mean winding angle as a 
function of I. Lower-right inset: semilogarithmic behavior of 
the winding. Dotted-dashed lines show the best fit according 
to Eq. (3), with « = 2.76 ± 0.1. The Solid line in the main 
frame shows the fit with n = 8/3 for comparision with SAWs. 
(b) The rescaled probability density function of the winding 
angle at lengths I = 50, 200, 400 in units of lattice spacing 
compared to the standard Gaussian density, 



of the random walk and of percolation clusters. In our 
case, since the dimension of the islands is 2, it suggests 
that they are compact unlike the clusters in critical per- 
colation. 

Statistical behavior of the coastlines here is similar to 
the statistics of rocky shorelines studied recently in [21]. 
The winding angle statistics of the shorelines is consis- 
tent with the prediction of SLE 8 / 3 . To be more serious 
about the similarities between the 2D-KPZ coastlines and 
rocky shorelines and, moreover, giving another justifica- 
tion for conformal invariancy, we compute the winding 
angle statistics for the contour lines of 2D-KPZ surface. 

Duplantier and Saleur studied in [22] the winding an- 
gle between the two endpoints of a finite SAW in two 



dimensions. Using Coulomb gas methods, they found 
that the distribution of winding angle is Gaussian with 
the winding variance of ~ (8/g) InL, where L is the dis- 
tance between the endpoints and g is Coulomb gas cou- 
pling parameter which is related to k by g = 4/k. They 
have also shown that the winding angle at a single end- 
point relative to the global average direction of the curve 
is a Gaussian with variance of (4/g)\nL. Wieland and 
Wilson found in [23] that the variance in the winding 
at typical points along the curve is 1/4 as large as the 
variance in the winding at the endpoints. 

We define the winding angle 9 (as used in [21]), as the 
angle between the line joining two points separated by a 
length I along the curve and the local tangent in the ref- 
erence point, measured counterclockwise in radian. For 
conformally invariant curves of diffusivity k, its variance 
behaves like [21] 



2 k 



■In/, 



(3) 



where the average is taken over an ensemble of 1800 
curves by moving the reference point along each curve. 

As shown in Fig. 5(a), the variance of winding angle 
for iso-height lines at a cut made at mean height has a 
logarithmic behavior. The best fit (dotted dashed lines), 
corresponds to k = 2. 76 ±0.1 which is compared with the 
fit by setting n = 8/3 for SAWs (solid line). The length 
scale , /, is measured in units of lattice spacing which is 
set to unity, on square sample size of 2048 2 . For each 
configuration the largest loop is selected, so within the 
length scale /, the curves do not have a preferred direc- 
tion and the mean winding angle is zero (upper-left inset 
in Fig. 5(a)). The rescaled probability density function 
of the winding angle at different lengths is shown in Fig. 
5(b), which converges to a standard Gaussian density. 
The results are consistent with [5]. We find no changes 
in the winding statistics of contour lines at different level 
sets. 

In conclusion, studying the statistics of iso-height lines of 
saturated 2D-KPZ surface at different level sets we find 
that the fractal dimensions of cluster heights and their 
perimeter remain unchanged when changing the height of 
the cut. We also find that the exponent associated with 
the distribution function of the cluster size (the mass of 
clusters is considered here) changes as a function of the 
height of the cut. It linearly decreases for cuts made 
below the mean height and crosses over to an almost lin- 
ear fluctuation around a specific value above the mean 
height. We also tested another exponents related to the 
distribution of length of the loops, the area of the loops 
and the radius of gyration which all change at different 
level cuts (the results are not included in this paper). 
The winding angle statistics of the contour lines also sug- 
gests that their statistics is comparable to that of SAWs. 
This confirms that the conformal invariant property of 
the contour lines are given by SLE 8 / 3 . 
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